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ABSTRACT 

The extinction-to-backscatter ratio (S a ) is an important parameter used in the determination of the aerosol extinction and 
subsequently the optical depth from lidar backscatter measurements. We outline the algorithm used to determine S a for 
the Cloud and Aerosol Lidar and Infrared Pathfinder Spacebome Observations (CALIPSO) lidar. S a for the CALIPSO 
lidar will either be selected from a look-up table or calculated using the lidar measurements depending on the 
characteristics of aerosol layer. Whenever suitable lofted layers are encountered, S a is computed directly from the 
integrated backscatter and transmittance. In all other cases, the CALIPSO observables: the depolarization ratio, 5, the 
layer integrated attenuated backscatter, p’, and the mean layer total attenuated color ratio, y, together with the surface 
type, are used to aid in aerosol typing. Once the type is identified, a look-up-table developed primarily from worldwide 
observations, is used to determine the S a value. The CALIPSO aerosol models include desert dust, biomass burning, 
background, polluted continental, polluted dust, and marine aerosols. 
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1. INTRODUCTION 

S a is an intensive aerosol property, i.e., a property that does not depend on the number density of the aerosol but rather 
on such physical and chemical properties as size distribution, shape and composition. These properties are governed 
primarily by the source of the aerosol. The accuracy of the S a value used in the lidar inversions is dependent in part, on 
the correct identification of the type of aerosol. The CALIPSO aerosol typing will utilize observation data (both 
climatological and field campaigns), knowledge of emission sources inferred from surfaces types, and CALIPSO 
aerosol measurements such as depolarization ratio and the ratio of the backscatter coefficients at 1064 nm and 532 nm, 
referred to as the color ratio, to identify aerosol type.. 

The importance of S a to lidar measurements is best illustrated by the solution of the single scattering lidar equation. 
While several solutions of the lidar equation are in use, we focus here on the Femald solution [1]. For illustrative 
purposes consider a ground based lidar system with a field of view so small that the multiple scattering effects can be 
neglected. The single scattering lidar equation for such a system can be written as, 

N( z) = T\—^— -yP T ( z)T 2 (z~zJAz 

hv 4n(z-z s ) , u 

where N(z) is the received signal photon count from a sample volume of thickness Az at altitude z, r\ is the system 
efficiency, At is the measurement integration period, P is the average laser output power, hn is the photon energy, h is 
Planck's constant, n is the optical frequency, A is the receiving telescope area, z s is the lidar site altitude, and (3 T (z) is the 
total volume backscatter coefficient. T 2 (z-zs) is the two-way atmospheric transmission over the range z-z s , 

T 2 (z — z s ) = exp [-2 \ z z a T ( l)dl] 

s 1 / 


where g t is the total atmospheric extinction coefficient. 

Both molecular and aerosol scattering contribute to the atmospheric extinction and backscatter coefficients, 
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cr T (z) = <J m (z) + <J a (z), 
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and 


P T (z) = P m (z) + P a (z) , j 4 

where the subscripts (m) and (a) refer respectively to molecular and aerosol scattering. Molecular scattering varies 
inversely to the fourth power of the wavelength (~X" 4 0117 ), whereas the wavelength dependence of aerosol scattering 
depends on the size distribution, shape, and refractive index of the aerosol particles. The total backscatter coefficient p T 
is calculated from the measured lidar profile using, 


X( z)exp 
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where X(z) = N(z)(z-z s ) 2 , N(z) is given by Eq. 1.1, and S a = a a /(3 a and S m = o m /p m are the extinction to backscatter ratios 
for aerosol and molecular scattering. S m is known from theory and is equal to 87i/3 sr. The reference altitude z 0 is 
chosen where the aerosol scattering is negligible, i.e. p a (z 0 )~0. The molecular backscatter coefficient profile is 
computed from theory using model atmosphere values for the temperature T(z) and pressure P(z). The value of S a on 
the other hand should be determined on a case by case basis and depends on the aerosol composition, size distribution, 
and shape. 

Several measurements have been made to determine S a for different aerosol types. These measurements, primarily at 
532 nm, were performed using a variety of techniques, including Raman lidar [2-4], slant path techniques, and 
transmittance methods [5], and a nephelometer modified to directly measure extinction and near 180-degree backscatter 
[6]. The assumption that S a has a fairly constant value throughout the vertical column was widely used in the lidar 
community for several years. Recent studies [7, 8] have shown that S a can vary significantly with altitude. 


2. LEVEL II EXTINCTION PRODUCTS 

CALIPSO will produce two level II aerosol extinction products depending on the choice of the aerosol extinction-to- 
backscatter ratio: an approximate product and a CALIPSO extinction product. The approximate extinction product will 
be derived from a fixed S a (35 and 30 sr at 532 and 1064 nm, respectively). These values are chosen because they 
correspond to S a of the clean rural or background aerosol. AERONET analyses described below show that clean 
background aerosol is a frequently encountered aerosol type in the atmosphere [9]. In addition, experience with LITE 
measurements shows that these values are not likely to cause the failure of the extinction calculation. The CALIPSO 
extinction product is derived from S a based on the best estimate developed using CALIPSO measurements and the most 
up-to-date field observations. This paper discusses methods of making this estimate and describes an algorithm to 
achieve this. 


3. S A FOR LOFTED AEROSOL LAYERS 

The transmittance method requires clear air above and below the layer so that the transmittance through the layer can be 
determined. Fig. 3.1 is an example of an aerosol layer lofted above a clear air region at 1km. The transmittance method 
uses the following equation describing the relationship between optical depth and integrated attenuated backscatter: 

r'=-^-{i-exp(- 

^ . 3.1 

Here y' is the integrated (from layer top to base) attenuated backscatter defined by, 
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/• base 

/= | PJr)Ta(r)dr 

J top 


t is optical depth, r| is a multiple scattering parameter. The quantities y\ and x describe characteristics of a feature, i.e., 
they are associated with the backscatter and/or extinction of particles only. Note that the effective two-way 
transmittance is T 2 = exp (-2r|x). If we define an effective S a , S* = r|S a , we can rewrite Eq. 3.3 as follows: 


2 / . 3.4 

The effective two-way transmittance is typically obtained by fitting the return both above and below a feature to a 
reference profile [5]. 



Fig. 3.1 The attenuated scattering ratios of a LITE return showing a case when the S a value can be calculated directly using the transmittance method 


4. S A SELECTION FROM LOOK-UP-TABLES 

One of the objectives of the algorithm is to estimate the appropriate value of S a within 30% of the true value. The 
strategy is to identify aerosol type and then use a look-up table to select values of S a and r| appropriate for the layer. The 
selection scheme uses the observed backscatter strength and depolarization to identify aerosol type, to the extent 
possible, from among one of the six types. In most cases, the depolarization is directly related to the hydration state of 
the aerosol. The backscatter and depolarization are not sufficient to fully constrain the model selection, however. 
Therefore, additional data is used to narrow down the choices of aerosol types based on the lidar observables. The 
selection algorithm uses the lidar observables and the International Geosphere-Biosphere Programme (IGBP) surface 
types. The land/water mask is based on the World Vector Shoreline (WVS) product for coastal information, and the 
Digital Chart of the World (DCW) for inland water areas. 


5. AEROSOL TYPES FROM AERONET 

Aerosol type is highly variable on time scales as short as a few hours [10]. Aerosol optical measurements must therefore 
be made at short time scales (about 30 minutes) to develop a large data base which can be used to derive statistically 
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significant correlations and from which type-specific characteristics can be deduced. The AERONET measurements are 
likely to provide such a data base albeit for total column rather than vertically-resolved measurements. Cluster analysis 
was used for grouping multi-year AERONET data set based on several optical and physical characteristics of the 
aerosol [9]. 

This study used whole AERONET archive (up to December 2002) of measurements and inversions to develop a type- 
specific set of mean optical properties of aerosols. Since the clustering algorithm we used requires a fixed number of 
clusters a priori , this number is determined by using successively higher numbers of clusters until no new significant 
clusters are formed. The data clusters formed in this way group all records that have statistically significant similarities 
in one category. Details of the analyses methods can be found in [9]. We discuss below the results of the cluster analysis 
for the aerosol types and the distributions of S a at 532 nm. The clustering analysis yields six distinct types of aerosols 
and each of these has distributions of the descriptive optical parameters used in the analysis. In this classification, 
polluted continental refers to a lightly-loaded soot-free pollution normally found in rural areas and is good 
approximation for rural aerosol. Dirty pollution refers to pollution containing significant amounts of absorbing species. 
Desert dust is assumed to be mostly mineral soil. Biomass burning is an aged smoke aerosol consisting primarily of soot 
and organic carbon. Polluted marine aerosol consists primarily of seasalt with traces of polluted continental species. 



Lidar Ratio (sr) 

Fig. 5.1 Distribution of all S a values at 532 nm 

Both urban and dirty pollution consist of the same species but the large imaginary part of the refractive index of dirty 
pollution suggests that this type of aerosol contains a significantly larger (than urban) fraction of soot. Fig 5.1 shows the 
distribution of S a for all records. Table I shows the mean and standard deviation (std. dev) of S a for the six aerosol types. 
An S a value of 32 ±6 sr for the clean northern hemisphere aerosol was measured during Aerosols99 [11]. A similar 
aerosol termed ‘background-like’ aerosol and originating on the European continent but devoid of any strong biomass 
or fossil fuel burning signature yielded an S a value of 35 sr [7]. The value estimated by this study for clean continental 
is slightly higher than these values mostly because the clean continental AERONET sites are not very pristine. 

The polluted marine S a value of 51sr is significantly larger than that measured during the SEAS experiment [12], and 
other estimates of S a for marine aerosols ranging between 20-30 sr [13-15]. 

The biomass burning value of 58 sr is not at variance with the measurements of Voss et al. (2001) of 60 ± 6 sr off the 
west coast of Africa, Ansmann et al. (2001) of 70 sr for biomass burning influenced aerosol advected from the Indian 
subcontinent during INDOEX. Recent studies of a limited climatology of 26 AERONET sites [16] at which biomass 
burning, coal combustion, urban/industrial, oceanic, and dust are the predominant types found S a values of 60±8 sr for 
biomass burning. 

In the case of desert dust, the AERONET retrievals used in the cluster analysis were based on the assumption of 
spherical particles. Dust particles are actually non-spherical, though, which has an unquantifiable effect on the accuracy 
of the results. The desert dust is lower than S a measurements by Voss et al. (2001) (41 ± 8 sr) using a Micropulse Lidar 
for African dust, Sasano and Browell (1989) (52 ± 10 sr) and measurements of Liu et al. (2002) of Asian dusts (42 - 55 
sr) found using a high- spectral-resolution lidar and a combined Raman elastic-backscatter lidar values. 

Measurements by Ansmann et al. (2001) at the Sagres island off the Portuguese coast showed the S a value for pollution 
emanating from continental Europe to vary between 50 and 70 sr. During INDOEX measurements, S a values of 
polluted continental aerosol originating from northern and northeastern part of India, known for high emissions of black 
carbon, were made by Franke et al. (2001). They found values ranging from 49 to 70 sr. Measurements of a stagnant 
airmass at Bondville (a polluted continental site) yielded S a values of 64 ± 4 sr [17]. 
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To check for the consistency of the individual categories, we divided the measurements in each category into five 
optical depth classes and plotted the S a of each class within a category. For each category, the magnitude of S a is fairly 
constant as shown in Fig. 5.2 below. The optical depth is an extensive property, i.e., a property that depends on the 
amount of aerosol. The Sa value on the other hand, is an intensive property and is relatively constant across optical 
depth classes within categories. This means that despite changes in the optical depth, the aerosol S a is consistently the 
same in each category. While this is not a sufficient validation of the clustering method, it is an indication of the ability 
of the algorithm to group similar data-sets using the prescribed variables. 

Distributions of the S a values by type (Fig. 5.3) show that these are mostly single mode Gaussian for polluted 
continental, biomass burning, polluted marine, and dirty pollution and somewhat bimodal for dust and clean continental. 
The S a values obtained from the AERONET measurements have several limitations including assumptions about the 
sphericity of the particles, homogeneity of the atmospheric column, and homogeneity of the aerosol layer. Despite these 
limitations, the study is useful in estimating (to first order) the variation of S a values and the corresponding phase 
functions by aerosol type and location. The phase functions of the six aerosol types are shown in Figure 4. This 
information is useful for the retrieval of optical depths from spaced-based lidar measurements such as the Cloud 
Aerosol Lidar and Infrared Pathfinder Spaceborne Observations (CALIPSO). 
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Fig. 5.2 Variation of the extinction/backscatter ratio over the optical depth within individual categories. For all the categories the variation with 
optical depth is small except for polluted continental and marine aerosols. 
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Fig. 5.3 Distributions of the extinction to backscatter ratios for the six aerosol types found by cluster analysis. 



Angle (deg.) Angle (deg.) Angle (deg.) 

Fig. 5.4 Phase functions of the six aerosol types 


Table 1 . Extinction to backscatter ratio statistics for the six aerosol types 


Biomass Clean Polluted Polluted Dirty 

Dust Burning Continental Continental Marine Pollution 

No. of Records 

22202 

26662 

20307 

55667 

6527 

12548 

Mean (sr) 

32 

58 

45 

58 

51 

83 

Std. Dev. (sr) 

14 

23 

20 

20 

19 

23 

Minimum (sr) 

10 

10 

10 

10 

10 

11 

Maximum (sr) 

146 

150 

149 

149 

147 

150 
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6. CALIPSO AEROSOL MODELS 






Fig. 6.1 The size distributions, and microphysical properties of the CALIPSO aerosol models. For each model, the S a (in sr) at 532 nm and 1064 nm, is 
shown in green and red numbers, respectively. r gjV and o g;V are the geometric mean radii and standard deviations of the distribution respectively. is 

the complex refractive index at wavelength i 


Some of the CALIPSO aerosol models were derived from AERONET clusters. Fig. 6.1 (a)-(f) show the physical and 
chemical properties of the CALIPSO aerosol models and the corresponding S a values at 532 nm (green) and 1064 nm 
(red). The AERONET cluster analysis yielded six distinct types of aerosol. Only three of these clusters (desert dust, 
biomass burning, polluted continental) were used to characterize the CALIPSO aerosol models. The CALIPSO model of 
background and marine aerosols are not derived from the AERONET measurements. 

The AERONET records of the background cluster had low mean optical depths (<0.05 at 673 nm). The microphysical 
properties derived from these are likely to have large uncertainties [18]. The CALIPSO background aerosol model (Fig. 
6.1c) was derived by fitting size distributions and refractive indices to measurements of S a of long-range continental 
transport [6]. Note that the S a values for this aerosol type are used to generate the approximate extinction product 
described in section 2 above. The AERONET marine aerosol cluster is comprised of a small number of records (< 4% 
of the total). The CALIPSO marine aerosol model (Fig. 6. If) is derived from the parameters measured during the SEAS 
experiment 12]. The CALIPSO polluted dust is a mixture of the AERONET desert dust and biomass burning clusters 
(Fig. 6.1e). 


7. S A SELECTION SCHEME 

The input parameters - the magnitude of attenuated backscatter, altitude, location, surface type, depolarization ratio, and 
mean attenuated backscatter coefficient measurements - are used to identify the type following one of eleven pathways 
in Fig. 7.1. In Fig. 7.1, pathway 1 is a lightly loaded aerosol layer found over snow/ice/tundra regions such as Antarctica 
and the clean Artie. Arctic haze, by virtue of the high integrated backscatter value, will be classified as polluted 
continental following pathway 2. Desert dust (Pathway 4) is expected to have a volume depolarization ratio greater than 
0.2 and should be the predominant selection pathway for layers with substantial fractions of non-spherical particles. 
Pathways 3 and 5 allow for mixing with biomass burning smoke which will depress the volume depolarization ratio to a 
value below 0.2. Pathway 6 is a clean non-desert land surface at which the aerosol loading is close to the background 
values. Pathway 7 is a highly polluted land surface such as would be found in urban areas. Pathway 8 accounts for 
continental pollution advected off the coast and entrained in the marine boundary layer. Pathway 9 is marine boundary 
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layer aerosol usually found in the deep ocean and consisting primarily of seasalt. Pathways 10 and 11 are elevated 
aerosol layers over land and ocean, respectively, of biomass burning smoke. Note that elevated dust layers are found in 
Pathway 4. 

Studies are underway to determine optimum threshold values of 5, (3’, and y to be used in the typing scheme. The values 
shown in Fig. 7.1 are initial estimates based on LITE measurements and, in the case of depolarization, on a limited set of 
observations and models [14, 19-22]. The goal is to base typing decisions on these observables as much as possible and 
avoid the use of geographic information. Therefore, the threshold values of 5, (3’, and y are implemented as runtime 
parameters that can be adjusted using a configuration script. When lofted layers are encountered under favorable 
conditions, S a is computed directly from the integrated backscatter and transmission. We expect the algorithm will 
evolve significantly and the CALIPSO measurements will provide a much larger set of 5, (3’, and y, measurements than 
is currently available. 


Launch Build Lidar Ratio Selection Flow Chart 



Fig. 7.1 Flowchart of the CALIPSO S a selection scheme for tropospheric aerosols. The values shown are the extinction-to -backscatter ratios at 532 nm 
and 1064 nm in parentheses. Note that the mean attenuated backscatter coefficient P' = y7z where z is the layer thickness 


8. CONCLUSION 

An algorithm for determining aerosol type from knowledge of emission sources inferred from surfaces types, and 
CALIPSO aerosol measurements (such as depolarization ratio and the mean attenuated backscatter coefficients at 532 
nm) has been developed. These aerosol types are characterized using observation data from both climatological and field 
campaigns. The strategy identifies aerosol type and uses a look-up-table to assign an extinction-to-backscatter ratio to 
the aerosol layer. 

A global data set, AERONET, has been used to identify main clusters of aerosol types and to determine microphysical 
properties of aerosol groups. This characterization is augmented by measurements where the uncertainty in the 
AERONET retrievals is high. New measurements from CALIPSO will significantly enhance this algorithm by making 
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available a vastly improved data base of 5, |3' , and x- Such data sets will be used to further refine the probability 

distribution functions and threshold values of 8, p' , and % used in the algorithms. The flowchart in Fig. 7.1 is therefore a 
preliminary selection scheme. This approach will ensure that the science mission requirements for optical depth accuracy 
of 40% error assuming a 30% uncertainty in the extinction-to-backscatter ratio will be realized. 
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